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Abstract. I discuss the latest progress in calculating baryon-baryon interactions with lattice QCD. The 
latest results from the NPLQCD collaboration for nucleon-nucleon scattering along with preliminary cal- 
culations of hyperon-nucleon scattering are presented. 
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1 Introduction 

Perhaps the greatest challenge facing those of us working 
in the area of strong interaction physics is to be able to 
rigorously compute the properties and interactions of nu- 
clei. The many decades of theoretical and experimental 
investigations in nuclear physics have, in many instances, 
provided a very precise phenomenology of the strong inter- 
actions in the non-perturbative regime. However, at this 
point in time we have little understanding of much of this 
phenomenology in terms of the underlying theory of the 
strong interactions. Quantum Chromo Dynamics (QCD). 

The ultimate goal is to be able to rigorously com- 
pute the properties and interactions of nuclei from QCD. 
This includes determining how the structure of nuclei de- 
pend upon the fundamental constants of nature. Any nu- 
clear observable is essentially a function of only five con- 
stants, the length scale of the strong interactions, ^qcd, 
the quark masses, m„, m^, nis, and the electromagnetic 
coupling, tte (at low energies the dependence upon the 
top, bottom and charm quarks masses is encapsulated in 
^qcd)- Perhaps as important, we would then be in the 
position to reliably compute quantities that cannot be ac- 
cessed, either directly or indirectly, by experiment. 

The only way to rigorously compute strong-interaction 
quantities in the nonperturbative regime is with lattice 
QCD. One starts with the QCD Lagrange density and per- 
forms a Monte-Carlo evaluation of Euclidean space Green 
functions directly from the path integral. To perform such 
an evaluation, space-time is latticized and computations 
are performed in a finite volume, at finite lattice spac- 
ing, and at this point in time, with quark masses that are 
larger than the physical quark masses. To compute any 
given quantity, contractions are performed in which the 
valence quarks that propagate on any given gauge-field 
configuration are "tied together" . For simple processes 
such as nucleon-nucleon scattering, such contractions do 
not require significant computer time compared with lat- 
tice or propagator generation. However, as one explores 



processes involving more hadrons, the number of contrac- 
tions grows rapidly (for a nucleus with atomic number A, 
charge Z and strangeness S*. the number of contractions 
is {A + Z)\{2A - Z- S)\S\ ED, and a direct lattice QCD 
calculation of the properties of a large nucleus is quite im- 
practical simply due to the computational time required. 
The way to proceed is to establish a small number of ef- 
fective theories, each of which have well-defined expansion 
parameters and can be shown to be the most general form 
consistent with the symmetries of QCD. Each theory must 
provide a complete description of nuclei over some range 
of atomic number. Calculations in two "adjacent" theories 
are performed for a range of atomic numbers for which 
both theories converge. One then matches coefhcients in 
one EFT to the calculations in the other EFT or to the 
lattice, and thereby one can make an indirect, but rigor- 
ous, connection between QCD and nuclei. It appears that 
four different matchings are required: 

1. Lattice QCD. Lattice QCD calculations of the prop- 
erties of the very lightest nuclei will be possible at some 
point in the not so distant future fT| . Calculations for 
A < 4 as a function of the light-quark masses, would 
uniquely define the interactions between nucleons up 
to and including the four-body operators. Depending 
on the desired precision, one could possibly imagine 
calculations up to A ~ 8. 

The chiral potentials and interactions have been de- 
termined out to the order where four-body interac- 
tions contribute [2] . As with the EFT constructions in 
the meson-sector and single-nucleon sector, the num- 
ber of counterterms proliferates with increasing order 
in the expansion, and at some order one looses rigorous 
predictive power without external input. Lattice QCD 
will make model independent determinations of these 
counterterms. 



^ Notice that the number of contractions in hypernuclei con- 
taining a small number of hyperons is less than in ordinary 
non-strange nuclei. 
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Exact Many-Body Methods. During the past decade 
one has seen remarkable progress in the calculation of 
nuclear properties using Green Function Monte-Carlo 
(GFMC) with the A^ig-potential (e.g. Ref. [3]) and 
also the No-Core Shell Model (NCSM) (e.g. Ref. g]) 
using chiral potentials. Starting with the chiral poten- 
tials, which are the most general interactions between 
nucleons consistent with QCD, one would calculate the 
properties of nuclei as a function of all the parame- 
ters in the chiral potentials out to some given order in 
the chiral expansion. A comparison between such cal- 
culations and lattice QCD calculations will determine 
these parameters to some level of precision. These pa- 
rameters can then be used in the calculation of nu- 
clear properties up to atomic numbers A ~ 20 — 30. 
The computer time for these many-body theories suf- 
fers from the same ~ (^!)^ blow-up that lattice QCD 
does, and for a sufficiently large nucleus, such calcula- 
tions become impractical. 

Another recent development that shows exceptional 
promise is the latticization of the chiral effective field 
theories [SJlEUail] • This should provide a model-independRltffi'^iP^'^*'^ ^^'^ meetmg 
calculation of nuclear processes once matched to lattice 
QCD calculations. 

Coupled Cluster Calculations. In order to move 
to larger nuclei, A< 100 a technique that has shown 
promise is to implement a coupled-clusters expansion 
(e.g. Ref. [10]). One uses the same chiral potential that 
will have been matched to lattice QCD calculations, 
and then performs a diagonalization of the nuclear 
Hamiltonian, after truncating the cluster expansion, 
which itself contains arbitrary coefficients. The results 
of these calculations will be matched to those of the 
NCSM or GFMC for A - 20 - 30 to determine the 
arbitrary coefffcients. This method is unlikely to be 
practical for very large atomic numbers. 
Very Large Nuclei and Density Functional The- 
ory (Maybe) To complete the periodic table one needs 
to have an effective theory that is valid for very large 
nuclei and nuclear matter. A candidate that has re- 
ceived recent attention is Density Function Theory (DFT) 
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Fig. 1. as a function of m^. This figure is reproduced with 
the permission of LHPC. 



mesons and the interactions between two mesons, primar- 
ily due to limitations in computational power, the last 
few years has seen an increasing number of precise calcu- 
lations of the properties of nucleons at the available pion 
masses. I wish to show two that will be of interest to the 



2.1 The Matrix Element of the Axial Current 

The matrix element of the axial current in the nucleon, 
gA, is a fundamental quantity in nuclear physics as it 
is related to the strength of the long-range part of the 
nucleon-nucleon interaction via PCAC. During the last 
year LHPC [T^ has for the first time calculated gA at 
small enough pion masses where Heavy Baryon Chiral Per- 
turbation Theory (HBxPT) should converge. The results 
of this calculation, previous calculations and the physical 
value, are shown in Fig.[T] Also shown is the chiral extrap- 
olation along with is uncertainty (the shaded region) . The 
physical value lies within the range predicted by chiral 
extrapolation of the lattice calculations. 



2.2 The ISIeutron-Proton Mass Difference 



(e.g. Refs. [TTlll2j ). It remains to be seen if this is in 
fact a viable candidate. There is reason to hope that 
this will be useful because there is clearly a density ex- 
pansion in large nuclei with a power-counting that is 
consistent with the Naive Dimensional Analysis (NDA) 
of Georgi and Manohar [T3] . The application of DFT 
to large nuclei is presently the least rigorously devel- 
oped component of this program. 
The latticized chiral theory mentioned previously can 
also be applied to the infinite nuclear matter problem. 
This work is still in the very earliest stages of explo- 
ration, but this looks promising |7j. 



During the last year it was realized that one could use 
isospin-symmetric lattices to compute isospin-breaking quan- 
tities to NLO in the chiral expansion [TS]. This is achieved 
by performing partially-quenched (unphysical) calculations 
of the nucleon mass [16 in which the valence quark masses 
differ from the sea-quark masses (those of the configu- 
rations) and determining counterterms in the partially- 
quenched chiral Lagrangian. The NPLQCD collaboration 
found that, in the absence of electromagnetic interactions, 
the neutron-proton mass difference at the physical value 
of the quark masses is 



M„ - Mp 



= +2.26 ± 0.57 ± 0.42 ± 0.10 MeV (1) 



2 Lattice QCD Calculations of Single 
Nucleon Properties 

While the lattice QCD calculations have historically fo- 



which is to be compared with estimates derived from the 



= -H2.05±0.3. 



P I physical 

We see that the lattice determination is consistent with 



Cottingham sum-rule of M„ 
We see that the lattice dete: 

what is found in nature. The lattice calculation is expected 
cused on the meson sector, both the properties of single to become significantly more precise during the next year. 
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Fig. 2. The tttt scattering length in the 1 — 2 channel as 
a function of m^/f^. The dashed curve corresponds to the 
unique prediction of tree-level xPT, while the shaded region 
is the fit to the resuhs of the NPLQCD calculations [21]. The 
NPLQCD calculations are performed at a single lattice spacing 
of 6 ~ 0.125 fm. 



3 Hadron Scattering from Lattice QCD 

To circumvent the Maiani- Testa theorem [T7] , which states 
that one cannot compute Green functions at infinite vol- 
ume on the lattice and recover S-matrix elements except at 
kinematic thresholds, one computes the energy-eigenstates 
of the two particle system at finite volume to extract the 
scattering amplitude [18] . The scattering amplitude of two 
pions in the / = 2-channel has been the process of choice 
to explore this technique, and to determine the reliability 
and systematics of lattice calculations. 



3.2 Ktt Scattering 

Studying the low-energy interactions between kaons and 
pions with K'^n~ bound-states allows for an explicit ex- 
ploration of the three- flavor structure of low-energy hadronic 
interactions, an aspect that is not directly probed in rnr 
scattering. Experiments have been proposed by the DIRAC 
collaboration 24J to study Kn atoms at CERN, J-PARC 
and GSl, the results of which would provide direct mea- 
surements or constraints on combinations of the scattering 
lengths. In the isospin limit, there are two isospin channels 
available to the Ktt system, / = ^ or / = | . The width of a 
K'^TT~ atom depends upon the difference between scatter- 
ing lengths in the two channels, F ~ (ai/2 "03/2)^7 (where 
ai/2 and 03/2 are the I = and I = ^ scattering lengths, 
respectively) while the shift of the ground-state depends 
upon a different combination, AEq ~ 2ai/2 + 03/2- 

The NPLQCD collaboration calculated the K^tt^ scat- 
tering length, af^+^+, in the same way as we computed 
TTTT scattering, requiring only the additional generation of 
strange quark valence propagators ^Sj. As ax+Tr+ was 
calculated (at a single lattice spacing of 6 ~ 0.125 fm) 
at three different pion masses, a chiral extrapolation was 
performed. This extrapolation depends upon two coun- 
terterms, one from the crossing-even and one from the 
crossing-odd amplitudes, and the important point is that 
their coefficients have different dependence upon the me- 
son masses. Therefore, both can be determined from the 
results of the lattice calculation, as shown in Fig. [3] and 
therefore, the scattering lengths in both isospin channels 
can be predicted, as shown in Fig. [3] This is a another 
demonstration of the combined power of lattice QCD and 
xPT. 



3.1 TTTT Scattering 

NPLQCD has calculated tttt scattering in the 1 — 2 chan- 
nel at relatively low pion masses using the mixed-action 
technique of LHPC. Domain- wall valence propagators are 
calculated on MILC configurations containing staggered 
sea-quarks. The lattice calculations were performed with 
the Chroma software suite [T^l^ on the high-performance 
computing systems at the Jefferson Laboratory ( JLab) . 

This calculation demonstrates the predictive capabili- 
ties of lattice QCD combined with low-energy EFT's. By 
writing the expansion of 771^02 as a function of m^//^, 
it has been shown that when inserting the values of 171^^ 
and /jr as calculated on the lattice. The lattice spacing 
effects in this mixed action calculation, which naively ap- 
pear at 0(6^), are further suppressed, appearing only at 
higher orders [22ll23j . The first non-trivial contributions 
from the partial-quenching in this calculation are shown 
to be numerically very small, and therefore a clean extrac- 
tion of the counterterm in the xPT Lagrangian is possible. 
The chiral extrapolation, as shown in Fig. [2l is found to 
have a smaller uncertainty at the physical point than the 
experimental value. 



3.3 Nucleon-Nucleon Scattering 

A few years ago we realized [^[?7] that even though the 
scattering lengths in the nucleon-nucleon system are un- 
naturally large, and much larger than the spatial dimen- 
sions of currently available lattices, rigorous calculations 
in the NN-sector could be performed today. There are two 
aspects to this. First, it is unlikely that the scattering 
lengths in the NN sector are unnaturally large when com- 
puted on lattices with the lightest pion masses that are 
presently available, m^^^ 250 MeV. Second, it is not the 
scattering length that dictates the lattice volumes that can 
be used in a Luscher-type analysis, but it is the range of 
the interaction, which is set by for the NN interaction. 
While the Luscher asymptotic formula are not applicable 
when the scattering length becomes comparable to the 
spatial dimensions, the complete relation is still applica- 
ble [H], 
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Fig. 3. The left panel shows the 95% confidence ellipses for the counterterms L5 and Lktv that contribute to Kn scattering at 
NLO in xPT as extracted from the mixed-action lattice QCD calculation of NPLQCD at a lattice spacing of 6 ~ 0.125 fm |25) . 
The right panel shows the 95% confidence ellipses for the Kn scattering lengths from a combination of a lattice QCD calculation 
and xPT Also shown are the 38% confidence ellipses from a Roy-Steiner analysis and from a 0{p'^) xPT analysis. 



where 6{p) is the elastic-scattering phase shift evaluated 

at kinetic energy T = 2 (^y/^PTJT^ - Mn^ , and the sum 

in eq. ^ is over all triplets of integers j such that | j | < yl 
and the limit yl ^ cx3 is implicit f27l . 

The same set of configurations that was used to pro- 
duce the scattering lengths for tttt and Kn scattering, dis- 
cussed above, were used to construct the NN correlation 
functions, and hence to extract the energy shifts between 
two nucleons at finite volume and twice the single nucleon 
mass at finite volume. At the pion masses accessible, the 
NN scattering lengths are found to be of natural size, set 
by the inverse pion mass. Only one of the pion masses 
is within the region described by the low-energy EFT's, 
and as such a prediction of the scattering length from lat- 
tice QCD alone cannot be made. However, when combined 
with the physical values of the scattering lengths, an al- 
lowed region for the scattering lengths as a function of 
pion mass can be made, as shown in Fig. [H 

3.4 Hyperon-Nucleon Scattering 

Computation of hyperon-nucleon (and hyperon-hyperon) 
scattering amplitudes with lattice QCD is performed in 
the same way as the computation of nucleon-nucleon scat- 
tering amplitudes. The required strange quark propaga- 
tor is not computationally expensive compared with the 
light-quark propagators, and so the calculation of these 
amplitudes costs very little more than that of the nucleon- 
nucleon amplitudes. NPLQCD has computed the correla- 
tion functions for hyperon-nucleon scattering in the chan- 
nels with interpolating fields AN (spin s = and s = 1), 
E~n (spin s = and s = 1), AA, S~ S" and 

is currently analyzing them in order to extract scatter- 
ing phase-shifts for - 290,350,500 and 600 MeV. In 
fig. [5] we show the preliminary correlation function and 
the effective scattering length for S^n in the spin-singlet 
channel for ^ 500 MeV and at a lattice spacing of 



h ~ 0.125 fm. Om preliminary analysis of these quantities 
gives a phase-shift of 5 — —61.8 ± 7.2° at an energy of 
T = 26.9 ± 7.2 MeV which would give a scattering length 
of a = +0.61 ± 0.12 fm for an energy-independent k cot 5. 

4 Hadronic Potentials 

Lattice QCD cannot, in general, isolate the potential be- 
tween hadrons, as the potential by itself is not an ob- 
servable quantity. An exception to this statement is when 
both hadrons are infinite massive, and consequently their 
spatial separation can be fixed (this of course generalizes 
to N-hadrons). Currently, Silas Beane, William Detmold, 
Kostas Orginos and I are performing a quenched calcula- 
tion of the potential between two B-mesons in the heavy 
quark limit. In this limit the light degrees of freedom (dof) 
of each B-meson decouple from the heavy-quark dynamics, 
and their spin, s;, becomes a good quantum number. The 
potential between the two B-mesons when constructed 
from xPT has the same form as that for NN, but with 
different values of the counterterms that enter. There- 
fore, calculating the potential between two B-mesons may 
provide qualitative insight into the NN-system. Our cal- 
culation is not the first attempt to extract the poten- 
tial between heavy hadrons, but is performed with a lat- 
tice spacing approximately half of that of the previous 
calculations [341135]. The quark mass was chosen to give 
~ 420 MeV and rup - 700 MeV. 

Fig. [6] shows our preliminary results for the central 
potentials as a function of B-meson separation (the ten- 
sor potential is found to be very small, and so we have 
also shown the potentials at r = \/2 and a/S which are 
not separated into tensor and central). One clearly sees 
short-distance repulsion in the channel with the quantum 
numbers of s-wave NN interactions, while one finds large 
attraction in the other channels. There is a correction due 
to the finite lattice spacing that is to be added to the re- 
sults of the lattice calculation, and the leading correction 
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Fig. 4. Scattering lengths in the and — ^Di NN channels as a function of pion mass. The experimental value of the 
scattering length and NDA have been used to constrain the extrapolation in both BBSvK [2811291150] and W [31II32II55] power- 
countings at NLO. 
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Fig. 5. The correlation function (left panel) and the effective scattering length plot (right panel) for neutron- X"" obtained in 
the coarse MILC lattices with L ~ 2.5 fm and ~ 500 MeV. 
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Fig. 6. Preliminary calculations of the potential between two B-mesons in the heavy quark limit. The vertical axis is the 
difference between the central potential at a separation r and that at r = 7b. The fine (red) dots correspond to the leading 
finite-lattice spacing correction that must be added to the potential obtained in the lattice calculation. 



is stiown in Fig. [S] Adding this factor gives the continuum 
potential between two infinitely massive B-mesons in the 
lattice volume, while the extrapolation to infinite volume 
potential has not yet been performed. The tensor inter- 
action between hadrons has not been directly isolated in 



previous works. In Fig.[7]we show our preliminary calcula- 
tion of the tensor potential between two infinitely massive 
B-mesons. We find a potential that is consistent with zero 
within the uncertainty of the calculation, indicating that 
the tensor potential in this system is Vr^ 40 MeV. 
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Fig. 7. Preliminary calculations of the tensor potential be- 
tween two B-mesons in the heavy quark limit. 



5 Outlook 

It is clear that lattice QCD is an important part of the 
future of nuclear physics, and perhaps is an even more im- 
portant part of the future of hypernuclear physics where 
there is significantly less experimental data to guide the- 
oretical constructions. Lattice QCD, when combined with 
effective field theory, is just starting to make rigorous pre- 
dictions of few-body observables, and we can expect sig- 
nificantly more progress as the computational resources 
dedicated to these calculations is increased. The formal 
tools are in place to explore the two-nucleon sector, and 
are being put in place to explore the hyperon-nucleon sec- 
tor [36 l l26 l [37] . all that is required is computational power. 
However, a coherent effort involving both numerical cal- 
culations and formal developments is presently required 
to move beyond the two-body sector and to explore hy- 
pernuclei. 

I would like to thank my collaborators in this work, 
Silas Beane, Paulo Bedaque, Kostas Orginos, William Det- 
mold, Tom Luu, Elisabetta Pallante and Assumpta Par- 
reno. 
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